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Aims. We present new radio continuum observations of ten BIMA SONG galaxies, taken at 1.4 GHz with the Very Large Array. These 
observations allow us to extend the study of the relationships between the radio continuum (RC) and CO emission to 22 CO luminous galaxies 
for which single dish CO images have been added to interferometric data. New Spitzer infrared (IR) images of six of these galaxies have been 
released. The analysis of these high resolution images allowed us to probe the RC-IR-CO correlations down to linear scales of a few hundred 
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ABSTRACT 



pc. 

O Methods. We compare the point-by-point RC, CO and mid-IR intensities across entire galaxy disks, producing radial profiles and spatially 
^ resolved images of the RC/CO and RC/mid-IR ratios. 

Results. For the 22 galaxies analysed, the RC-CO correlation on scales from ~ 10 kpc down to ~ 100 pc is nearly linear and has a scatter 
of a factor of two, i.e. comparable to that of the global correlations. There is no evidence for any severe degradation of the scatter below the 
^ , kpc scale. This also applies to the six galaxies for which high-resolution mid-IR data are available. In the case of NGC 5194, we find that the 

' non-thermal radio spectral index is correlated with the RC/FIR ratio. 
^ ' Conclusions. The scatter of the point-by-point correlations does not increase significantly with spatial resolution. We thus conclude that we 
^ 1 have not yet probed the physical scales at which the correlations break down. However, we observe local deviations from the correlations in 
- - ' regions with a high star formation rate, such as the spiral arms, where we observe a flat radio spectrum and a low RC/FIR ratio. In the intra-arm 
regions and in the peripheral regions of the disk, the RC/FIR is generally higher and it is characterized by a steepening of the radio spectrum. 
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1. Introduction 

One of the major goals in studying spiral galaxies is to un- 
derstand the relationship between the star formation process 
and the physical conditions in the interstellar medium. Since 
the discovery that stars form in molecular clouds, several ef- 
forts have been directed towards the study of the relation be- 
tween the emission of the CO molecule, which traces the bulk 
of molecular gas, and the other star formation indicators. 

Global studies have revealed relationships between the ra- 
dio continuum (RC), far-infrared (FIR) and CO emissions. The 
global correlation between the FIR and centimeter-wavelength 
RC emission is one of the strongest coiTelations in extragalac- 
tic research. Early works on this topic were reported by Helou 
et al. (1985) and de Jong et al. (1985). Condon (1992) and 
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more recently Yun et al. (2001) found that these two emis- 
sions are linearly coiTelated over five orders of magnitude in 
luminosity, with an rms scatter of less than a factor of 2. On 
global scales the FIR emission is well correlated with the CO 
emission (e.g. Devereux & Young 1990) and the CO emission 
is well correlated with the RC (Rickard et al. 1977; Israel & 
Rowan-Robinson 1984; Adler et al. 1991). 

What is most extraordinary about the FIR-RC or CO-RC 
correlations is that they couple emissions arising from com- 
pletely different processes. The usual explanation of the FIR- 
RC correlation invokes massive-star formation which accounts 
for the thermal radio emission via ionizing stars, for the non- 
thermal radio luminosity via supernova events, and for the 
FIR luminosity by means of massive stars heating the dust 
(Wunderlich & Klein 1988). However, this basic scenario has 
too many steps and too many parameters to explain the tight- 
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ness of the observed correlations. More complex models have 
been put forward, see Murgia et al. (2005, hereafter Paper I) for 
a brief summary, but a 'definitive' model has not yet emerged. 

Although there is still no consensus on the causes of the 
RC-FTR-CO correlations, nor any adequate explanation for 
their tightness, these correlations have been used to determine 
empirically the star formation rate (SFR) (Cram et al. 1998) 
and the star formation efficiency (SFE) in galaxies (Gao et al. 
2001, 2003; Murgia et al. 2002; Rownd & Young 1999). More 
generally it is found that these correlations are well described 
by a "Schmidt Law" of the type SFR oc p^^, where pgas is the 
total gas density, and the exponent A'^ typiciiUy ranges from 1.3 
to 1.5 (Kennicutt 1998). 

Detailed studies of individual objects at high spatial resolu- 
tion are essential to understanding whether there is a common 
physical process behind the observed correlations and possibly 
synthesizing the various empirical descriptions. 

An improved approach to the study of the spatial properties 
of the FIR-RC correlation was presented in Marsh & Helou 
(1995). They studied the FTR-RC correlation on intermediate 
scales of ~ 1-3 kpc and they found a small but systematic 
decrease in the FIR-RC ratio as a function of increasing ra- 
dial distance from the center. In contrast, Murgia et al. (2002) 
found that the ratio between the RC and CO emission is con- 
stant, to within a factor of 3, both inside the same galaxy and 
from galaxy to galaxy, down to kpc size scales. Leroy et al. 
(2005) extended this result to a sample of dwarf galaxies, find- 
ing that this relationship between CO and RC is the same in 
all galaxies. In this context, dwarf galaxies appear to be simple 
low-mass versions of large spirals. 

So far the low resolution of FIR observations has prevented 
the comparison of the FIR and RC emission on kpc scales for 
all but a few Local Group sources. Xu et al. (1992), studying 
the Large Magellanic Cloud found evidence for a breakdown 
in the FIR-RC correlation at a hnear scale of ~ 70 pc. In the 
Milky Way, Boulanger & Perault (1988) analysed the Orion 
region and showed that the FIR-RC correlation breaks down 
on the scale of a few hundred parsecs around a region of star 
formation. 

The situation has recently changed with the advent of the 
satelhte Spitzer that provides the necessary resolution of IR 
observations to extend the study of the correlation to sub- 
kpc scales in many external galaxies. These observations can 
shed new light on the study of the correlation on local scales. 
Recently, Murphy et al. (2006) presented an initial look at the 
far-infrared-radio correlation within the star-forming disks of 
four nearby galaxies using Spitzer images. 

An ideal sample to investigate the spatially resolved 
FIR-RC-CO correlations is the BIMA SONG sample. The 
Berkeley-Illinois-Maryland Association Survey of Nearby 
Galaxies (BIMA SONG) is an imaging survey of the CO (J=l- 
0) emission in 44 nearby spiral galaxies with a typical angular 
resolution of 6" (Regan et al. 2001 ; Heifer et al. 2003). For half 
of the galaxies in the BIMA sample single-dish CO images at 



55" resolution obtained with the NRAO 12 m telescope' have 
been added to interferometric data. The combined images rep- 
resent a significant improvement with respect to the previous 
CO extragalactic surveys, since in each galaxy a large area can 
be mapped at high resolutions avoiding the missing flux prob- 
lem. 

Recently we extended the study of CO-RC correlation 
down to a spatial resolution of ~ 100 pc, by collecting 20 cm D 
and B-array data of a subsample of nine bright BIMA SONG 
galaxies available in the VLA archive (Paper I). For these nine 
galaxies it has been possible to compare CO and RC images at 
two different angular resolutions: a low resolution (~ 55") set 
obtained from the NRAO 12m and VLA D-array, and a high 
resolution (~ 6") set. The point-by-point RC and CO bright- 
nesses across the entire galaxy disks have been compared and 
the ratio of the CO and RC flux has been defined as a way to 
characterize the CO-RC correlation. 

In this work we present the study of 13 more BIMA SONG 
galaxies which allow us to complete the study, started with 
Paper I, of 22 CO-luminous BIMA galaxies. 

In Section 2 we summarize the observations used in this 
paper and the data reduction procedure. Section 3 describes 
briefly how we performed the data analysis. The results are pre- 
sented in Section 4 for the whole sample, and are discussed in 
Section 5. We analysed the 24 /zm-RC correlation for six galax- 
ies of the sample, for which high resolution Spitzer mid-IR im- 
ages are available, at a spatial resolution ranging from ~400 pc 
to ~100 pc. We present the results of this analysis in Section 
6. In Section 7 we compare the expectations of the leaky box 
diffusion model model, proposed in Paper I, with the observed 
trend of spectral index versus RC-FTR ratio for NGC 5194. In 
Section 8 we summarize our results. 

2. Observations and data reduction 

We present new radio observations of ten BIMA SONG galax- 
ies, taken at 1.4 GHz with the Very Large Array synthesis tele- 
scope in its B configuration. These observations allow us to 
extend the study of the relationships between the radio contin- 
uum and CO emission, started in Paper I, to 22 CO-luminous 
galaxies, for which single dish CO images have been added to 
interferometric data. 

To these newly observed we add three more BIMA SONG 
galaxies: NGC 3351, NGC 4303 and NGC 5457 whose 
archival data were available, and NGC 1068, for which we use 
the high resolution image of the FIRST public survey of the 
VLA. The galaxies studied in this work are Usted in Table 1, 
along with their main geometrical parameters. A summary of 
the VLA-B array data, including date and length of observa- 
tions, and beam size and rms level of produced images for each 
target is reported in Table 2. 

We selected a bandwidth of 50 MHz for the two IF chan- 
nels centered at 1.465 and 1.385 GHz. The second of the two 
IFs was corrupted by interferences and had to be discarded, re- 

' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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Source 



a(/2000) (5(72000) 



PA. RC3 Type, Nucleus 



name 




(h m s) 




C ") 


(Mpc) 


n 


(°) 




NGC 


0628 


01 36 41 


.7 


+15 46 59 


7.3 


24 


25 


SA(s)c 


NGC 


1068 


02 42 40, 


.7 


- 00 00 48 


14.4 


33 


13 


(R)SA(rs)b, Syl.8 


NGC 


2903 


09 32 10 


.1 


+21 30 02 


6.3 


61 


17 


SAB(rs)bc,Hn 


NGC 


3351 


10 43 31, 


.1 


+ 11 42 14 


7.4 


40 


13 


SB(r)b,Hn 


NGC 


3521 


11 05 48, 


.0 


- 00 02 04 


7.2 


58 


164 


SAB(rs)bc,HlI/L2 


NGC 


3627 


11 20 15, 


.1 


+12 59 22 


11.1 


63 


176 


SAB(s)b, T2/Sy2 


NGC 


3938 


11 52 49, 


.6 


+44 07 14 


17 


24 





SA(s)c,Hn 


NGC 


4303 


12 21 54, 


.9 


+04 28 25 


15.2 


27 





SAB(rs)bc, T2 


NGC 


4321 


12 22 54, 


.8 


+ 15 49 20 


16.1 


32 


154 


SAB(s)bc,T2 


NGC 


4826 


12 56 43, 


.6 


+21 40 59 


4.1 


54 


111 


(R)SA(rs)ab,T2 


NGC 


5248 


13 37 32, 


.06 


+08 53 07 


22.7 


43 


110 


SAB(rs)bc, HH 


NGC 


5457 


14 03 12, 


.48 


+54 20 55 


7.4 


27 


40 


SAB(rs)cd 


NGC 


7331 


22 37 04, 


.09 


+34 24 56 


15.1 


62 


172 


SA(s)b,T2 



Coordinates, distances, inclinations, position angles, RC3 type and nuclear classifications from Heifer et al. (2003). Inclination is defined 
assuming i=0° for face-on galaxies. Position angles are measured in degrees from north through east. For nuclear spectral classification: H 
n = H II nucleus, Sy = Seyfert nucleus, L = LINER, and T=transition object, which has [O I] strengths intermediate between those of H II 
nuclei and LINERs 



Table 2. VLA B-array observation summary 



Source 


Date 


Duration 


Beam 




name 




hours 


HPBW(") 


(juJy beam"') 


NGC 0628 


03-NOV.-2003 


1.5 


5.7x5.2 


40 


NGC 2903 


04,06-Nov-2003 


1.0 


5.9x5.3 


35 


NGC 3351 


28-Jul-1986' 


1.8 


4.5 X 4.1 


46 


NGC 3521 


04,06-Nov-2003 


1.0 


4.9 X 4.8 


65 


NGC 3627 


04,06-Nov-2003 


1.0 


5.8x5.7 


35 


NGC 3938 


04,06, lO-Nov-2003 


2.0 


5.8x4.9 


50 


NGC 4303 


13-Sep-1982** 


1.0 


6.5 X 4.6 


100 


NGC 4321 


04,06,10-Nov-2003 


1.2 


5.8x5.1 


30 


NGC 4826 


06,10-Nov-2003 


1.0 


5.5 X 5.2 


45 


NGC 5248 


04,06, lO-Nov-2003 


1.8 


6.1 X 5.2 


50 


NGC 5457 


06-Aug-1990*" 


1.0 


5.8x5.3 


50 


NGC 7331 


03-NOV.-2003 


1.0 


5.3 X 4.8 


50 



a) Final noise of total intensity images. 



' Archival data: program A073; 
" Archival data: program COND; 
•** Archival data: program AVI 80; 



ducing the signal to noise ratio by a factor of V2. Nevertheless, 
the time spent on each source was sufficient to guarantee that 
the dynamic range of RC images is larger than that of the cor- 
responding CO images. 

The visibility data were calibrated, Fourier transformed, 
and deconvolved with the NRAO package AIPS following stan- 



dard procedures. Self-calibration was applied to minimize the 
effects of phase variations. The final cleaned maps were con- 
volved with the beams listed in Table 2. The flux densities were 
brought to the scale of Baars et al. (1977) using the flux cali- 
brators 3C48, 3C138 and 3C286. 
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Fig. 1. (a) Radio contours of galaxies overlaid on the optical images (grey scale), taken from the Digitized Sky Surveys (DSS) archive. See 
Table 2 for sensitivity and beam sizes of radio images. Contours start at 2 cr, and increase by steps of a factor of 2. 



In Fig. 1 we present the high-resolution radio image over- 
laid on the optical image for each galaxy, except for NGC 1068, 
for which we use the FIRST image. We can recognize differ- 
ent morphologies of radio emission, which are described in the 
notes on individual galaxies reported in Section 3.1. 



To improve the sensitivity of the radio continuum images to 
large-scale structure, D and C array data were combined with 
the B array data. When u-v data were available in the VLA 
archive we combined them in the visibility plane; otherwise D 
images were added in the image plane using the linear combi- 
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Fig.l. (i)) See Fig. l.a. 



nation method implemented in the MIRIAD task IMMERGE 
(Sault et al. 1995). 

We convolved the original RC and CO maps to common 

beam sizes for all but a few cases in which the data sets had 
intrinsic beam areas matching within 10%. As in Paper I, we 
have CO and RC images at two different angular scales: low 



resolution ( 55") and high resolution ( 6"). The image param- 
eters are listed in Tables 3 and 4, respectively. The beam sizes 
are listed in angular and Unear scales. 
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Fig. 2. (a) NGC 0628. Radio continuum and CO comparison in BIMA SONG galaxies. (Upper left) RC contours overlaid on CO emission 
for 55" images. {Lower left) RC contours overlaid on CO emission for 6" images. Contours start at 5cr and scale by a factor of V2. (Upper 
right) Ij 4 as a function of Ico for the 55" and 6" data sets. The insets show the box grids used; the grid size corresponds to the resolution, 
so that each point is essentially independent. The dotted (dashed) lines show the mean (dispersion) of the Ij 4 - Icq relation obtained for the 
sample containing 22 galaxies (this work and Paper I). The solid blue line is a weighted fit to the points which takes into account the errors in 
both coordinates. (Lower right) Radial profiles of Ico and I1.4 (top) and Qco/rc = Ico/Ii.4 (bottom) for the 55" and 6" images. The mean and 
dispersion of qco/Rc = log(2co/Rc) is also given. 



3. Data analysis 

Following the data analysis performed in Paper I, we compared 
the point-by-point RC and CO brightnesses across the entire 
galaxy disks. We overlaid regular grids of rectangular beam- 
sized boxes on both the RC and CO images, omitting areas with 



obvious background continuum sources. We averaged all pixel 
values within each box to calculate the 1 .4 GHz RC brightness 
(I1.4 ) and the CO integrated intensity (Ico )■ 

Figure 2 shows the results obtained for the sample of thir- 
teen galaxies. For each galaxy, the upper and bottom left pan- 
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Galaxy 


VLA data 




O-RC 






*) 






("X") 


(mjy bm"') 


("X") 


(kpcxkpc) 


(Jy bm"' km s"') 


NGC 0628 


D 


60x60 


0.30 


55x55 


2.0x2.0 


10.1 


NGC 1068 


D 


60x60 


0.20 


55x55 


3.8x3.8 


31.2 


NGC 2903 


D 


60x60 


0.20 


55x55 


1.7x1.7 


36.1 


NGC 3351 


D 


54x54 


0.30 


55x55 


1.9x1.9 


33.6 


NGC 3521 


D 


60x60 


0.15 


55x55 


1.9x1.9 


44.9 


NGC 3627 


D 


54x54 


0.20 


55x55 


2.9x2.9 


55.1 


NGC 3938 


D 


54x54 


0.15 


55x55 


4.5x4.5 


18.4 


NGC 4303 


D 


54x54 


0.40 


55x55 


4.0x4.0 


19.6 


NGC 4321 


D 


54x54 


0.30 


55x55 


4.3x4.3 


37.5 


NGC 4826 


D 


60x60 


0.10 


55x55 


1.1x1.1 


26.5 


NGC 5248 


D 


54x54 


0.15 


55x55 


6.0x6.0 


28.3 


NGC 5457 


D 


48x48 


0.25 


55x55 


2.0x2.0 


24.7 


NGC 7331 


D 


48x48 


0.25 


55x55 


4.0x4.0 


47.0 



a) The Unear resolution refers to the CO beam. 

b) cTco is rms level in undipped image of integrated intensity. 



Table 4. High resolution CO and RC image parameters. 



Galaxy 


VLA data 




O-RC 






6) 






("X") 


(mJy bm"') 


("X") 


(pcxpc) 


(Jy bm~' km s 


NGC 0628 


B+D 


7.2x5.3 


0.035 


7.2x5.3 


250x190 


1.3 


NGC 1068 


B 


8.8x5.5 


0.3 


8.8 x5.5 


614x384 


5.3 


NGC 2903 


B+D 


6.8x5.7 


0.03 


6.8x5.5 


210x170 


4.6 


NGC 3351 


B+D 


7.3x5.2 


0.05 


7.4x5.2 


265x187 


3.6 


NGC 3521 


B+D 


8.8x5.7 


0.05 


8.8x5.7 


310x200 


5.5 


NGC 3627 


B+D 


7.3x5.7 


0.03 


7.3x5.8 


390x310 


2.8 


NGC 3938 


B+D 


5.8x4.9 


0.03 


5.9x5.4 


480x440 


1.9 


NGC 4303 


B+D 


7.3x5.4 


0.1 


7.3x5.4 


538x398 


2.1 


NGC 4321 


B+D 


7.2x5.3 


0.03 


7.2x4.9 


562x382 


2.5 


NGC 4826 


B+D 


7.5x5.4 


0.04 


7.5x5.4 


150x100 


4.9 


NGC 5248 


B+D 


6.9x5.7 


0.06 


6.9x5.7 


759x627 


2.1 


NGC 5457 


B+D 


5.8x5.3 


0.04 


5.7x5.4 


210x190 


1.6 


NGC 7331 


B+C 


6.1x4.9 


0.03 


6.1x4.9 


440x360 


3.8 



a) The linear resolution refers to the CO beam. 

b) (Tco is rms level in undipped image of integrated intensity. 



els show the overlay of 1 1.4 (contours) on Ico (color) images, 
respectively, at low and high resolution. 

The upper right panels plot I1.4 versus Ico at low (soUd 
dots) and high (open dots) resolution. In the bottom right cor- 
ner insets, the upper and lower panels show the box grids for 
the low and high resolution images, respectively. Only points 
above 2cr have been plotted. 

The bottom-right panel shows the radial profiles of Ii 4 and 
Icq (top) and their ratio Qco/rc = Ico(Jy sr ' km s"')/ I1.4 
(mJy sr~^ ) (bottom). The value of the mean and dispersion of 
qco/RC = log(2co/Rc) is given on the bottom right plot of the 
radial variation in Qco/rc- For further details on the analysis 
technique see Paper I. 



3.1. Comments on individual galaxies. 

NGC 0628. This galaxy, seen almost face-on, has a clear grand- 
design spiral structure. The CO emission is strongly confined 
to the spiral arms (Regan et al. 2001). Our RC high-resolution 
image (Fig. 1) is deficient in radio emission, except for a bright 
spot corresponding to an HE region (Kennicutt & Hodge 1980) 
and an X-ray source (Colbert et al. 2004). The low-resolution 
RC image shows emission extended over about 6'. The point- 
by-point correlation for this source is flat and the radial profile 
shows a decrease of the qco/RC ratio. 

NGC 1068. This is a nearby and luminous Seyfert 2 galaxy, 
which has been extensively studied from radio wavelengths to 
X-rays. The nuclear emission is clearly dominated by the radio 
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Fig. 2. (b) NGC 1068. See Fig. 2a for technical details. 



source related to the AGN. For the study of the point-by-point 
correlation we do not use the low resolution data because the 
source is unresolved. We show only the point-by-point corre- 
lation and the radial profiles measured in the 6" resolution im- 
ages. Recently, Spinoglio et al. (2005) found that the infrared 
emission is due to both an AGN component and to the starburst 
component in the circum-nuclear ring of ~ 3 kpc in size. In the 
BIMA image the CO strongest emission is not associated with 
the nucleus of the galaxy but has a ring-like distribution around 
it. We omitted values from the fit of the correlation 1 1.4 and Ico 
measured in a central region of ~ 15" where the AGN emis- 
sions of the compact core and prominent lobes dominate (Ho 



& Ulvestad 2001). The resulting correlation is steeper than the 
mean of all galaxies. 

NGC 2903. This large starburst galaxy has a large-scale bar 
visible in CO and in high-resolution RC images. The Ha im- 
ages (Knapen et al. 2002) show star-formation along the bar 
and in the arms, mostly in the southern arm, and in the circum- 
nuclear region. The bar has been studied in detail by Sheth 
et al. (2002), who have compared the distribution of the molec- 
ular gas, traced using CO emission, and star formation activity, 
traced using the Ha emission line. They have found that the CO 
emission is brightest along the leading edges of the bar, while 
the majority of the Ha emission is offset, with values between 
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Fig. 2. (c) NGC 2903. See Fig. la for technical details. 



200-800 pc towards the leading side of the CO. In the point-by- 
point RC-CO correlation (upper-right panel of Fig. 2 (c)) it is 
possible to distinguish between data points of the central 30", 
which show a steeper correlation, and those of the bar emis- 
sion, with a flatter correlation. The same trend is not evident 
in the radial plot (lower-left panel) where the ratio is averaged 
over annuli in which the bar-like structure is invisible. 

NGC 3351. This is a typical ringed galaxy. Conselice et al. 
(2000) suggest that the nuclear starburst in this galaxy is due 
to the dynamical influence of the prominent bar. In the BIMA 
image the gas appears concentrated nearly entirely along the 
"twin peaked" area of the nuclear region (Kenney et al. 1992). 
Because of its small size (< 1'), the galaxy is not completely 



resolved in the low resolution image. The few low-resolution 
data are aligned with the extrapolation of the high-resolution 
data to lower intensity. 

NGC 3521. This highly tilted Sbc galaxy shows strong and 
broad Ho- absorption around the nucleus. This absorption com- 
ponent is most likely caused by a stellar absorption line of the 
central bulge component, characteristic of A-type stars. In or- 
der to have such a bright A-type star cluster, a large-scale star- 
burst may have occurred during the last 10^ yr (Sofue et al. 
1998). The high-resolution CO image shows a pronounced lack 
of emission at the center. Our high-resolution radio image (Fig. 
1) shows the spiral galaxy with its flocculent structure, and a 
strong background radio source to the southwest. The presence 
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Fig. 2. (d) NGC 3351. See Fig. 2a for technical details. 



of this strong background source degrades the dynamic range 
of the image and prevents us from reaching a better signal-to- 
noise ratio in the image. 

NGC 3627. This is an interacting spiral galaxy with a weak 
bar A warp is present in the optical image and is probably 
caused by tidal interaction with the other member of the Leo 
Triplet (NGC 3628). X-ray observations with Beppo SAX and 
Chandra suggest that the bulk of X-ray emission has a star- 
forming origin (Georgantopoulos et al. 2002). The molecular 
gas is observed in a central concentration along the leading 
edges of the bar, at the end of the bar, and extending over nearly 
4' in spiral arms that extend from the ends of the bar (Regan 
et al. 2001). The same unusual distribution is present in the RC 



emission at high resolution (see Fig. 1). The qco/RC radial pro- 
file (see lower right panel of Fig. 2f) shows the presence of a 
hook in the central region, corresponding to the nucleus, which 
exhibits a composite LINER/H II spectrum. Filho et al. (2000) 
reported the presence of a variable, flat-spectrum radio source. 

NGC 3938. For this galaxy the RC contours in the bottom- 
right panel of Fig. 2 (g) start from 3 cr, since the RC emission at 
high resolution is extremely weak, with steps of a factor of V2. 
It is a nearly face-on late-type spiral galaxy with well defined 
arms. It produces very faint radio emission, mostly distributed 
along the spiral arms (Fig.l). The CO emission is also faint, so 
studying the correlation is hard, given the low dynamic range 
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Fig. 2. (e) NGC 3521. See Fig. la for technical details. 



of both emissions, nevertheless the mean value is not far from 
the mean of the overall sample. 

NGC 4303. This is a nearby double-barred galaxy, which 
harbors a low-luminosity AGN. The molecular gas properties 
in the central kiloparsec have been studied by Schinnerer et al. 
(2002), who have found that the bulk of the molecular line 
emission comes from two straight gas lanes that run north- 
south along the leading sides of the large-scale primary bar 
Because of the low signal-to-noise ratio of our high resolution 
RC image, obtained from archive data, we can see an analogu- 
ous distribution of the radio emission. The RC-CO point-by- 
point correlation is linear, and the radial plot of qco/RC shows 



the presence of the central hook, corresponding to the nuclear 
region, classified as a Seyfert 2/LINER . 

NGC 4321. This is a well-studied, grand-design spiral 
galaxy, located within the Virgo Cluster. A barred radio struc- 
ture is present in our radio image (Fig. 1). The same distribution 
is visible in the high-resolution CO image. This is evident also 
from the point-by-point correlation and from the radial plot. 

NGC 4826. The 'Evil Eye' galaxy is a relatively isolated 
Sab spiral with a conspicuous dust lane. Ho et al. (1997), 
analysing spectroscopic parameters, have classified the nucleus 
as a transition object, which has a radio counterpart in our 20 
cm image (Fig.l). Archival observations at higher resolution, 
performed with VLA-A array, revealed the presence of clumps 
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Fig. 2. (f) NGC 3627. See Fig. la for technical details. 



of emissions. In the upper- and lower-right panels of Fig. 2 (1) 
we can see that the point-by-point and also the radial profile 
show a good correlation between the RC and CO emission. 

NGC 5248. This is a nearby barred galaxy with a spectac- 
ular grand-design spiral structure. Our radio image shows the 
presence of a barred structure (Fig. 1). The properties of the 
circum-nuclear region and the relation to the host galaxy have 
been discussed by Knapen et al. (2002). Massive star forma- 
tion, as traced by Hor emission, is very strong in the circum- 
nuclear region and in the (non-active) nucleus, and is located 
mainly along the spiral arms and in the bar. The RC-CO cor- 
relation seems to lie in the mean of the sample, and the radial 
profile does not show any significant trends. 



NGC 5457. This barred spiral galaxy is nearly face-on and 
tidally interacting with NGC 5474, NGC 5477 and Holmberg 
IV, which results in a distortion of the HI distribution in its 
outer regions (Davies et al. 1980). In our map (see Fig. 1) it 
presents a lopsided RC emission. Due to its large size (~ 10'), 
we could infer a loss of flux from extended regions (>3') us- 
ing the B-array, but the same asymmetric distribution of radio 
emission is evident in the D-array image (see the top right panel 
of Fig. 2 (n)). Kuntz et al. (2003) have studied the diff^use X- 
ray emission, which traces the spiral arms and is roughly cor- 
related with the Ha and far-UV emission. The X-ray source 
population has been studied with ROSAT HRI and PSPC ob- 
servations by Wang et al. (1999), who have suggested that the 
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Fig. 2. (g) NGC 3938. See Fig. 2a for technical details. 

detected sources result from recent massive star formation in 
the outer spiral arms. For this unusual source the point-by-point 
correlation is flatter than the mean of the sample. 

NGC 7331. This is a relatively highly inclined (/ = 77°) 
spiral galaxy with a ring visible in CO (e.g. Regan et al. 
2001), radio continuum (Cowan et al. 1994), submillimeter 
(Bianchi et al. 1998) and mid-infrared images (Smith 1998). 
Our high-resolution image clearly shows the ring-like structure 
and the central emission due to a very low-luminosity radio 
core (Cowan et al. 1994), which make this galaxy a prototype 
of a submillijansky AGN (Filho et al. 2004). The RC-CO cor- 
relation has a very low dispersion and the qco/RC is almost 
constant within the ring structure, which is about 1 kpc thick. 



4. Results 

We start by summarizing the results of the study of RC-CO 
correlation from kpc to sub-kpc scale for the 13 BIMA SONG 
galaxies analysed in this work. We will then consider the whole 
sample composed of the nine galaxies from Paper I and these 
13, in order to improve the statistical significance of our results. 

4.1. Point-by-point and radial correlations of the RC 
and CO emissions 

The point-by-point comparison of the 1 1.4 - Ico correlations are 
shown at both low (> 1 kpc) and high (< 1 kpc) resolution for 
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Fig. 2. (h) NGC 4303. See Fig. 2a for technical details. 



each galaxy in the upper-right panels of Fig. 2. The values of 
the slope and of the rms, obtained by a least squares fit to the 
measured values, are also given. 

We find that for four galaxies (namely NGC 0628, NGC 
1068, NGC 3938 and NGC 5457) the correlation is far from 
perfect. In the case of NGC 0628, the low-resolution images 
show that the RC emission is more extended than the CO one, 
at high resolution RC emission is faint and diffuse, while CO 
is tightly confined to the spiral arms. This drop of RC emis- 
sion causes a deviation from the mean Ii 4 - Ico relation, result- 
ing in a very flat slope. NGC 1068 is a very luminous Seyfert 
galaxy. Only the high-resolution correlation is studied here, 
since in low-resolution images the source is unresolved. The 



emission from the central region of ~ 15" is dominated by the 
compact core and prominent lobes associated with the AGN 
(Ho & Ulvestad 2001). Omitting these points from the fit (see 
dashed box in Fig. 2b), the remaining points have a larger dis- 
persion, and the resulting correlation is steeper than the mean. 
For NGC 3938, due to the very low dynamic range of both the 
RC and CO emission, all points are clustered in the Ii 4 - Ico 
plot. NGC 5457 is particularly peculiar since it exhibits a lop- 
sided RC emission, which is not seen in the CO high-resolution 
image. This galaxy is ~ 10' in diameter, so D-array observa- 
tions do not suffer from losses of flux emitted from extended 
regions, while in B-array observations emissions from regions 
more extended than 3' are lost. Nevertheless, both low and 



R. Paladino et al.: Thermal and non-thermal components of the interstellar medium at sub-kiloparsec scales in galaxies 



15 







g 




3a 3a Aa 


is 





Hi ■ 



1 «■ 



■\aa a . -a- 



JO'" 



ADD 9~ A flS^ 



■ ■ ■ I ■ 



^ 4 B 



Fig. 2. (i) NGC 4321. See Fig. la for technical details. 



high resolution images show an asymmetric distribution of ra- 
dio emission, probably due to the tidal interactions with other 
galaxies. The slope and dispersion of the point-by-point corre- 
lation are higher than the mean. 

The remaining nine galaxies (see Section 3.1 for individual 
features of the correlation of each galaxy) show slopes and dis- 
persions of the point-by-point correlation consistent with the 
mean values for the overall sample. 

In the bottom-right panels of Fig. 2 the radial profiles of 
Ii.4, Icq and of their ratio are shown. The mean and dispersion 
in qco/RC = log(Ico/Ii.4) are given in the bottom-right corner. 
Omitting NGC 1068 from the calculation, the average qco/RC 
for the remaining twelve sources is 1.20+ 0.25; the disper- 



sion in qco/RC within the individual sources ranges from 0.08 
to 0.49. The findings of the point-by-point correlation are re- 
flected in the radial analysis: the four anomalous galaxies show 
a large deviation from the mean behaviour in the radial profile. 
In particular, NGC 0628 shows a radial decrease of qco/RC of 
about one order of magnitude from central to external region. 
In contrast, NGC 1068 shows the opposite trend, but this is due 
to the fact that in the radial analysis the central region has not 
been omitted. 

Figure 3 shows the images of the ratio qco/RC for each 
source. All images are shown on a common grey-scale, so that 
color variations from source to source indicate real variations in 
qco/RC- The qco/RC images of NGC 0628, NGC 3938 and NGC 
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Fig. 2. G) NGC 4826. See Fig. la for technical details. 

5457 clearly show their deficiency of RC emission; NGC 5457 
is particularly remarkable with its lopsided RC emission de- 
spite the almost symmetric CO distribution. In the qco/RC im- 
age of NGC 1068 the enhanced RC emission from the central 
AGN is evident. There are some differences in the distribution 
of qco/RC among the other sources. For NGC 3351, NGC 5248 
and NGC 7331 there is very little variation in qco/RC- This is 
also evident from their low dispersions (Fig. 2). In NGC 3521, 
NGC 4321 and NGC 4826, variations that follow the galaxies' 
spiral structure are present. The ratios are enhanced along the 
dust lanes in the bars of NGC 2903 and NGC 3627. 



4.2. Results for the overall sample 



In Fig. 4 we present histograms of the RC-CO correlations for 
the sample of 22 BIMA SONG galaxies. The left panel shows 
the histogram of the RC-CO correlation slopes. Although the 
slope varies from one galaxy to another (in the range of 0.29 
to 1.6, including also two peculiar galaxies: NGC 0628 and 
NGC 1068) the mean value is 1.1± 0.3. In the right panel we 
show the histogram of qco/RC obtained from the point-by-point 
measures of Icq and I; 4. Averaged over all points of the grids, 
<qco/Rc> = 1.2 + 0.3, i.e. the uncertainty in the mean CO/RC 
ratio across all galaxies is less than a factor of 2. 
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Fig. 2. (m) NGC 5248. See Fig. 2a for technical details. 



In Fig. 5 we show the correlation between RC brightness 
and CO integrated intensity of the 22 galaxies, corrected for 
the inclination. The points correspond to values computed as 
annular averages, i.e. they are radial profiles. We omitted all 
points of NGC 1068 dominated by the AGN radio emission 
from the fit. The resulting correlation has a slope of 1.1 and 
a rms of 0.3 dex, i.e. the RC-CO correlation is Unear with a 
dispersion of less than a factor of two. The dispersion increases 
progressively below an RC brightness of 10** mJy/sterad. In this 
range there are three anomalous galaxies (NGC 0628, NGC 
3938 and NGC 5457) : as we have noted in previous sections, 
their CO emission is clearly visible, while the RC emission at 
1 .4 GHz is so weak that the spiral structure is not detected at 



all at high resolution (see Fig. 1). These results confirm the 
tightness of the CO-RC correlation found in Paper 1, although 
there can be significant variations from source to source. 



The main result of this work is that for the sample of 22 
BIMA SONG galaxies we find the point-by-point and radial 
correlations between the CO and RC emissions to be consis- 
tent with the results of Paper 1. The CO-RC correlation extends 
over four orders of magnitude, and has a dispersion of less than 
a factor of 2. This appears to be a general characteristic for 
spiral galaxies whose radio emission is dominated by starburst 
activity. 
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Fig. 2. (n) NGC 5457. See Fig. 2a for technical details. 



5. Discussion 

5.1. The RC-CO decorrelation scale(s) 

One of the main aims of this work was to find out whether 
the RC-CO-FIR correlation persists at sub-kpc scales in galax- 
ies. The spatial resolution of our images is adequate for this 
purpose. In fact, we can explore the RC-CO correlation from 
scales of 10 kpc down to 100 pc. In order to estimate quanti- 
tatively the behaviour of this correlation we consider the rms 
scatter and the Pearson coefficient of correlations for the 22 
galaxies of the sample. We know that the relativistic electrons 
diffuse from their birth places. Therefore one expects that the 
spatial correlation of the emissions should break down below 



the characteristic diffusion scale-length of the radiating elec- 
trons, Ldiff ~ Vdift- ■ T, where v^is is the diffusion velocity and t 
is the electron lifetime. If we use the Alfven velocity (typically 
~ 10 km s"') in the warm ISM phase of the galactic disk, which 
may be responsible for turbulent diffusion of cosmic rays, as 
Vdiff, we get Ldiif « 1 kpc. This value is extremely uncertain, 
depending on the diffusion mechanism of electrons, which is 
not fully understood and can be considered a lower limit to the 
scale-length (if we assume as the typical wind velocity in star- 
burst galaxies v^ ~ 2000 km s ' as the diffusion speed, we 
obtain a much higher value of Ldiff). 

Figure 6 shows the Pearson coefficient and the scatter of 
the correlation versus the linear scale of the observations. For 
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Fig. 2. (o) NGC 7331. See Fig. 2a for technical details. 



each galaxy we plot two points, corresponding to the high and 
low resolution correlations, respectively. The median of the rms 
scatter increases slightly from 0.15 in the kpc range to 0.25 in 
the sub-kpc range, while the median value of the Pearson cor- 
relation coefficient (r) shows the same trend, decreasing from 
0.93 to 0.84. Thus, this further quantitative estimate of the 
goodness of the correlation confirms the result found consider- 
ing the scatter of the correlation, i.e. we cannot recognize any 
abrupt variation that could indicate a possible spatial scale of 
the decorrelation. 

In Fig. 6 we have labelled points corresponding to a 
Pearson coefficient smaller than 0.7 and to a scatter higher than 
0.4. We recognize those galaxies with an anomalous behaviour 



as has been emphasized in Section 4 and in Paper I: IC 342, 
NGC 0628, NGC 1068, NGC 4258, NGC 4736 and NGC 5457. 

The little degradation of the correlation at high resolution 
could be due to physical effects, such as the electrons' diffu- 
sion from their location of acceleration, or to observational ef- 
fects (some galaxies are not completely resolved by the low- 
resolution observations, and this can artificially increase the 
tightness of their correlation on kpc scales). 

We emphasize that our study of a statistically significant 
sample of 22 BIMA SONG galaxies does not reveal any ir- 
refutable decorrelation of the RC and CO emission on sub-kpc 
scales. This result indicates either that we have not yet probed 
the spatial scales at which the correlations break down or that 
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Fig. 4. Le/f Histogram of measured slopes for the overall sample of galaxies. Right Histogram of all measured qco/RC grid boxes for the 22 
sample galaxies. The means and dispersions are also given. 

5.2. Radial behaviour 
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Fig. 6. Rms scatter (top panel) and Pearson coefficient (bottom panel) 
of the RC/CO correlation for 2 1 galaxies versus linear scale of obser- 
vations. Points corresponding to a Pearson coefficient smaller than 0.7 
and to a scatter higher than 0.4 are labelled. 



there is a mechanism of regulation that compensates for the 
electrons diffusion. 



As noted in Paper 1 the most significant deviation from the 
close correlation between the FIR and RC emission is a small 
and generally monotonic decrease in FIR/RC as a function 
of radius. Marsh & Helou (1995), with a study at a resolu- 
tion of found a decrease of the FIR-RC ratio with the 
radial distance from the centers of galaxies, which they in- 
terpreted as a smearing of the RC due to the propagation of 
energetic electrons. The study of Marsh & Helou has seven 
galaxies in common with our sample: NGC 0628, NGC 2903, 
NGC 4303, NGC 5055, NGC 5194, NGC 6946 and NGC 7331. 
These galaxies show monotonic decreases in FlR/RC over 10- 
20 kpc in radius. In contrast, with the sample of galaxies in this 
work we study the radial distribution of the lco/Ii.4 ratio in the 
inner 10 kpc of the disks and we do not see any systematic trend 
with radius (see bottom right panels of Fig. 2). NGC 0628 is an 
exception: as we can see in Fig. 2 (a) there is a radial decrease 
of qco/RC analogous to that observed in the FlR-RC ratio by 
Marsh & Helou (1995). 

For the remaining galaxies, we have seen that the most co- 
herent structures in qco/RC are found along the spiral arms, with 
an enhanced qco/RC along the spiral arms and lower values in 
the interarm regions. 

Recent Spitzer FIR imaging of four nearby galaxies re- 
ported by Murphy et al. (2006) show that the spiral structure 
is visible in infrared/radio ratio images, which show enhanced 
values of FIR emission along the arms, with local peaks cen- 
tered on H 11 regions and depressed ratios located in the quies- 
cent interarm and peripheral regions of each galaxy. Gordon 
et al. (2004) find very similar features in Spitzer images of 
M 81. These images confirm the prediction of Paper 1, that 
when averaged in azimuthal rings, the fraction of the total area 
in an annulus subtended by the spiral arms is much larger at 
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Fig. 5. Correlation between the RC and CO emission in 22 galaxies. I1.4 and Icq have been corrected for the galaxy inclinations. Each point 
represents the average value within an annulus. The solid line is a weighted fit to the points shown, which takes into account the errors in both 
coordinates and has a slope of 1.1. The points of NGC 1068 clearly due to the AGN radio emission are omitted from this plot. Different symbols 
represent the galaxies studied in this paper, crosses are the galaxies of the sample investigated in Paper I. 



small radii than at large ones. It may be that the measured ra- 
dial falloff in the FIR/RC ratio in the IRAS and ISO observa- 
tions is primarily a consequence of the low resolution of the 
FIR observations. 

6. 24//m-RC high-resolution correlation with Spitzer 
data 

With the launch of the Spitzer Space Telescope (Werner et al. 
2004), it is now possible to map many large galaxies in the far- 
IR with unprecedented spatial resolution and sensitivity. The 
MIPS (Multiband Imaging Photometer for Spitzer; Rieke et al. 
2004) observations at 24 //m, 70 //m and 160 pm, have PSFs 
with angular resolutions of ~5.1", ~16"and ~38", respectively. 

In this section we focus on six galaxies in our sample for 
which Spitzer images have been recently released: NGC 3351, 
NGC 3521, NGC 3627, NGC 4826, NGC 5 1 94 and NGC 7331. 
With 6" angular resolution we reach spatial scales ranging from 
~400pcto~ 120 pc. 

In Fig.s 7 and 8 we show the point-by-point correlations be- 
tween the RC, CO and mid-IR 24 pm emissions for these galax- 
ies, sorted by increasing linear resolution. From these plots we 
cannot recognize any systematic variation of correlations from 
lower to higher linear resolution. The scatter and slope of these 



correlations are represented graphically in Fig. 9 as functions 
of the spatial resolution. There is no significant variation of the 
slope and the scatter of the correlations at this spatial resolu- 
tion. All three correlations are comparably tight: they extend 
over two orders of magnitude and have a scatter of less than a 
factor of two. What is clearly seen is that the 24/im-RC corre- 
lation has the smaller scatter and somewhat flatter slope, while 
the scatters of RC-CO and 24//m-CO correlations are compa- 
rable (see Table 5). 

In conclusion, it seems that the 24jum-RC-CO correlations 
persist at spatial scales of the order of ~ 100 pc. The observed 
small scatter in galaxies for which we obtain the higher spa- 
tial resolution indicates that we have not yet probed the spatial 
scales at which the correlations break down. An increase of the 
scatter is observed for the galaxy NGC 3627, which is char- 
acterized by a complex morphology. Probably the presence of 
bar and spiral arms increases the scatter of the high-resolution 
correlations. 

7. NGC 5194: a test for the leaky box model 

NGC 5194 has been recently observed with Spitzer. Calzetti 
et al. (2005) found that the 24 pm luminosity appears to be 
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Fig. 7. RC-CO, RC-24//m and 24ytim-CO correlations of NGC 7331, NGC 3627 and NGC 352 1 at 6" resolution. The linear resolution achieved 
for each galaxy is indicated in the lower-right of the right-most panels. 



closely correlated to sites of active star formation; an ini- 
tial look at the 24 jum-RC correlation has been presented by 
Murphy et al. (2006). We use these new FIR high-resolution 
images to test the predictions of the model proposed in Paper I. 
We have proposed that the CO-RC-FIR correlations can be ex- 
plained as a consequence of regulation by hydrostatic pressure. 
In high-pressure regions, a higher fraction of the gas becomes 
molecular, which enhances the CO emission. In these regions 
of strong FIR emission the higher pressure is also manifested 
in a stronger magnetic field, which enhances the synchrotron 
emission. 

In this model, the ratio between the confinement time, tc, 
and the radiative lifetime, fsyn, of the synchrotron electrons is 
related to the ratio between RC and FIR emissions (see Eq. 8 in 
Paper I). In particular, we expect that an increase of the RC/FIR 
ratio will result in a steepening of the synchrotron spectrum. 



Using the 70 jum Spitzer image of NGC 5 194 it is now pos- 
sible to study the point-by-point relation between the spectral 
index and the RC-FIR ratio. We measured the radio continuum 
spectral index between 1 .4 and 4.9 GHz using the images by 
Klein et al. (1984). These images have a resolution of ~ 16", 
comparable to that of the 70-yum observations, and incorporate 
short-spacing single-dish observations. The relation between 
the radio continuum spectral index a, defined as ly ~ v"", 
and the RC/FIR ratio is shown in Fig. 10. As the RC/FIR ratio 
increases by an order of magnitude, a monotonic steepening 
of the spectral index is observed. This behaviour is expected 
assuming a leaky-box model of cosmic rays in which the in- 
jection rate of the relativistic electrons. A, is proportional to 
the FIR and inversely proportional to the electron confinement 
time. The line shown in Fig. 10 represents the fit of the spectral 
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Fig. 8. RC-CO, RC-24/im and 24;um-CO correlations of NGC 3351, NGC 5194 and NGC 4826 at 6" resolution. 



index trend expected for a leaky box model with an injection 
spectral index of ao=0.5. 

The exact spectral shape has to be computed numerically. 
However, we have two asymptotic power law regimes; for tc <s; 
fsyn the synchrotron emissivity is: 

iv = /(«o)^^M) ■ B'+'^V-"" (1) 

while for tc » fsyn the spectrum is a power law with a spectral 
index a - ao + 0.5. 

The constants C/ and Cy depend only on fundamental phys- 
ical constants: C/ = 2.3 x 10"^^ erg Gauss"', Cy = 6.3 x 10'^ 
erg"^ Gauss"' s"' (see e. g. Pacholczyk 1970). 

The term /(a) has to be calculated numerically and for 
a'o=0.5 is /(a) ^ 3, while No-A-tc- Here A^o is the normaliza- 
tion of a cosmic -ray electron energy spectrum N(e) - N[)e^^, 
where N(e)de is the number density of cosmic -ray electrons 



with energies between e and e + de. For ao^O.S, No and A are 
expressed in erg ■ cm""*, and in erg ■ cm"^ s"', respectively. 

We tentatively relate the rate of energy injection in cos- 
mic ray electrons A to the FIR emissivity, and hence to the 
SFR, introducing a proportionality constant rf, so that: 

A^r] ;fir (2) 

The electron synchrotron lifetime at a given frequency is, in 
c.g.s. units: 

f.v« = 0.93-%^^-B-'V-5 

Substituting all of this into Eq. 1 we have the RC and FIR emis- 
sivity related by: 

- 14.. jj jjao-0.5.^-ao-0.5 Jc_ j-^^ 

7fIR tsyn 
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Galaxy 


Res 


Slope 


Scatter 


Slope 


Scatter 


Slope 


Scatter 




(pc) 


RC-CO 


RC-24;um 


24;um-CO 


NGC7331 


400 


0.72 


0.12 


0.63 


0.08 


0.83 


0.12 


NGC 3627 


350 


1.18 


0.27 


0.64 


0.16 


0.86 


0.26 


NGC 3521 


300 


0.53 


0.1 


0.47 


0.07 


0.71 


0.15 


NGC 3351 


230 


1.34 


0.2 


0.84 


0.16 


1.21 


0.26 


NGC 5194 


200 


1.28 


0.26 


0.56 


0.17 


0.77 


0.26 


NGC 4826 


120 


1.46 


0.21 


1.14 


0.13 


1.03 


0.16 


Average 




1.1±0.2 


0.19±0.02 


0.7+0.1 


0.13±0.01 


0.9±0.1 


0.20±0.02 



Table 5. Parameters of the RC-CO, RC-24yum and 24^m-CO correlations at 6". The scatter is expressed in dex. 
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Fig. 9. Rms scatter (top panel) and slope (bottom panel) of RC-CO, 
RC-24//m and 24yum-CO correlations versus linear resolution for six 
galaxies of our sample for which Spitzer images have been released. 
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Fig. 10. Point-by-point comparison between the spectral index and 
RC-FIR ratio, measured at a resolution of ~ 16". Points above 3 cr 
have been plotted. The line represents the fit of the leaky box model 
(see text). 



For our observations at v=1.4 GHz, assuming an injection 
spectral index of ao=0.5, we obtain: 

V IpiR lobs y'piR ' syn 

Expressing Eq. 2 in terms of FIR brightness, the value of tj 
allows us to infer the average electron injection rate at a given 
region of the disk: 

4.1x10-3*^ Ifir _i _3 

A - rj ■ erg s cm (6) 

Idisk /kpc MJy sterad ' 

where /disk is the thickness of the disk, and IfjR is the FIR 
brightness of the galaxy in units of MJy sterad 

In the case of NGC 5194, from the fit shown in Fig. 10, 
we estimate the value of ?/ - lO^s"'. We measure an av- 



erage FIR brightness of 25 MJy sterad which results in 

10"^^ erg s^'cm"^. 

If the CR electrons are injected in the energy range from 
€i and 62, this requires a injection power of Wf/=A In ^, for 
ao = 0.5. If we consider an electron energy spectrum ranging 
from y 1 to 10^, where y is the electron's Lorentz factor, it 
results in a power W^/ ^ 10"^^ erg s^'cm"^. 

This energetic requirement can be compared to the SFR, 
assuming that a fraction ^ of the total energy of a supernova 
(SN) goes into CR electrons, providing a yield of ^ ■ 10^' erg 
per SN. In a volume of 1 kpc-* a SN rate of v^ai^I SN yr"' 
yields 

WsFR - f 1 .4 ■ 10-^' y,'" , erg s"' cm-^ (7) 

vr A knr ^ 
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Condon (1992) calibrated the relation between the star forma- 
tion and the SN rate: 



( ysN \ 
\yr-i j 



0.041 



SFR(M>5Mq) 



Mgyr-i 



Hence we obtain: 



WsFR=^5.7-10 



-23 



SFR(M>5Mq) 



. Moyr"ikpc~3 



erg s cm 



(8) 



.(9) 



The star formation rate of NGC 5194 estimated from the 

extinction-corrected Ha luminosity is SFR ~ 4.2 M© yr~' 
(Scoville et al. 2001). Assuming a diameter of 10 kpc and a 
thickness of the disk of 1 kpc gives Wsm = ^ ■ 7.2 x 10"^^ 
erg s"' cm"-'. This value is comparable to our estimate if a 
fraction f ^ 10 of SN energy goes into the acceleration of 
CR electrons. 

The diffusive shock acceleration mechanism in supernova 
remnants is beheved to produce the majority of Galactic cos- 
mic rays (see Hillas 2005, for a recent review). Nevertheless, 
the understanding of the electron injection is poor. Typically 
it is assumed that 10% of the SN explosion energy goes into 
CR acceleration and of this 1% into the electron component. 
Recently Dyer et al. (2001) studying TeV emission found that 
the current energy content in relativistic electrons in SN 1006 
is about 7 xlO''^ ergs, i.e. 10"^ of the SN explosion energy. On 
the other hand, for Tycho's SNR Voelk et al. (2005) found a 
proton injection rate of 10 which, assuming an electron to 
proton ratio of ~ 0.01, means an electron acceleration rate of 
10"^. Thus our estimate for the conversion factor f is plausible 
from the point of view of the energy budget given the observed 
results on individual SNRs. 

8. Summary and conclusions 

1. We obtained new 1.4 GHz radio continuum images at 6" 
resolution for a sample of 13 nearby spiral galaxies belong- 
ing to the BIMA SONG survey. 

2. We analysed the behaviour of the point-by-point correla- 
tion between the RC and CO intensities from kpc to sub- 
kpc scales in these objects, extending the previous study by 
Murgia et al. (2005) to 22 galaxies. 

3. For most galaxies we find that the RC and CO intensities 
are nearly Unearly correlated (slope ^ 1.1 ± 0.3) over four 
orders of magnitude, with a scatter of a factor of 2 (0.3 dex). 
However, for four galaxies, namely NGC 0628, NGC 1068, 
NGC 3938 and NGC 5457, the correlation is not seen. NGC 
1068 is a Seyfert galaxy whose RC emission is strongly 
contaminated by the AGN. NGC 0628, NGC 3938 and 
NGC 5457 have low CO and RC surface brightness. The 
deviation from the RC-CO correlation is related to the fact 
that the RC emission is fainter than average at 6" resolu- 
tion. 

4. In the inner 10 kpc of the disks of the 22 galaxies analysed 
in this work the qco/RC ratio is constant with radius. The 
few exceptions are again NGC 0628, NGC 1068 and NGC 
3938. 

5. In order to determine whether the goodness of the RC- 
CO correlation degrades going from kpc to sub-kpc scales. 



we studied the variation of the correlation scatter and the 
Pearson coefficient as a function of the spatial resolution of 
our images. We find that the median scatter increases from 
0.15 to 0.25, going from ~10 kpc to 0.1 kpc. A similar trend 
characterizes the correlation coelTicient from about 0.93 to 
0.84. We conclude that the point-by-point RC-CO correla- 
tion deteriorates slightly when increasing the spatial reso- 
lution. However, we caimot identify any characteristic scale 
of decorrelation. 

6. We present the point-by-point 24 yum-RC and 24 //m-CO 
correlations at sub-kpc scales for six galaxies of our sam- 
ple. We probe the correlations on spatial scales ranging 
from ~ 450 pc to ~ 120 pc. The RC-CO and 24jum-CO cor- 
relations are comparably tight, with a scatter of ~ 0.2 dex. 
The 24yum-RC correlation has a scatter of only 0.1 dex and 
is somewhat flatter than the FTR-CO and RC-CO correla- 
tions. Moreover, we find that the scatter of the correlations 
does not increase with increasing spatial resolution. 

7. For the galaxy NGC 5194, using the Spitzer 70 //m image, 
it has been possible to test the leaky box diffusion model 
proposed in Paper I. A point-by-point comparison between 
spectral index and RC-FTR ratio, measured at a resolution 
of ~ 16", shows that regions with a flat spectral index have 
a reduced RC-FIR ratio. 

Acknowledgements. We thank the referee U. Klein for useful com- 
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